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Abstract. We have obtained discrete Eourier power spectra of a sample of persistent low-mass neutron-star X-ray 
binaries using long-term light curves from the All Sky Monitor on board the Rossi X-ray Timing Explorer. Our 
aim is to investigate their aperiodic variability at frequencies in the range 1 x 10“^ — 5 x 10“® Hz and compare 
their properties with those of the black-hole source Cyg X-1. We find that the classihcation scheme that divides 
LMXBs into Z and atoll sources blurs at very low frequencies. Based on the long-term (~ years) pattern of 
variability and the results of power-law fits {P{y) oc v~°‘) to the 1 x 10“^ — 5 x 10“® Hz power density spectra, 
low-mass neutron-star binaries fall into three categories. Type I includes all Z sources, except Cyg X-2, and the 
atoll sources GX9-I-1 and GX13-I-1. They show relatively flat power spectra (a < 0.9) and low variability (rms< 
20%). Type H systems comprise 4U 1636-53, 4U 1735-44 and GX3+1. They are more variable (20% < rms"^ 
30%) and display steeper power spectra (0.9 < a 1.2) than Type I sources. Type HI systems are the most 
variable {rms > 30%) and exhibit the steepest power spectra (a > 1.2). The sources 4U 1705-44, GX354-0 and 
4U 1820-30 belong to this group. GX9+9 and Cyg X-2 appear as intermediate systems in between Type I and II 
and Type H and III sources, respectively. We speculate that the differences in these systems may be caused by 
the presence of different types of mass-donor companions. Other factors, like the size of the accretion disc and/or 
the presence of weak magnetic helds, are also expected to affect their low-frequency X-ray aperiodic varibility. 
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1. Introduction 

A low-mass X-ray binary (LMXB) contains a neutron star 
which is accreting material via Roche lobe overflow from 
a companion star of spectral type later than A. Due to 
the high angular momentum of the accretion flow an ac¬ 
cretion disc is formed around the compact object (see e.g. 
White 1989). Unlike high-mass X-ray binaries, whose opti¬ 
cal spectrum is dominated by the emission from the mas¬ 
sive companion, in LMXBs the mass-losing star is usu¬ 
ally not seen owing to the contribution of the accretion 
disc (Cowly et al. 1991, Shahbaz et al. 1996). Lewin & 
van Paradijs (1985) proposed that the different classes of 
LMXBs might reflect the type of companion, with the 
brighter systems (usually those located in the galactic 
bulge) containing an evolved star and the fainter systems 
a main-sequence star. Based on infrared spectroscopic ob¬ 
servations, Bandyopadhyay et al. (1999) suggested that 
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not only the luminosity class but also the spectral type 
may be different. 

Various schemes have been proposed to categorise the 
LMXBs: X-ray spectral behaviour as a function of inten¬ 
sity (Parsignault & Grindlay 1978), cluster analysis of a 
large number of source characteristics (Ponman 1982), 
detailed X-ray spectral hts (White & Mason 1985), X- 
ray hardness-intensity and colour-colour diagrams (Schulz 
et al. 1988), age and location in the Galaxy (Naylor & 
Podsiadlowski 1993). In general, a bimodal distribution of 
sources and a correlation with luminosity is found. 

More relevant for the purpose of this paper is the clas¬ 
sification scheme in terms of the rapid aperiodic vari¬ 
ability and the patterns that these sources display in 
X-ray colour-colour diagrams (Basinger & van der Klis 
1989). In this scheme LMXBs are divided into two differ¬ 
ent subclasses, known as Z and atoll sources. In the fre¬ 
quency range 10“^-10^ Hz the power spectra of atoll and 
Z sources are represented by a power law with index 1- 
1.5 (atoll sources) and 1.5-2 (Z sources), which describes 
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Table 1. List of sources 


Source 

number 

Source 

name 

Alternative 

name 

Glass'" 

Lx xl0^“ 

erg s“^ 

r 

AP 

q3 

xVdoP 

rms" 

^orb 

hr 

Type I 

1 

4U 1617-15 

Sco X-1 

Z,-,I,2,B 

0.7 

920 

2.1 

0.67+0.06 

23/24 

14.4+0.2 

18.9 

2 

4U 1642-45 

GX 340+0 

Z,-,II,-,B 

2.2 

30 

2.0 

0.63+0.05 

43/24 

13.8+0.2 


3 

4U 1702-36 

GX 349+2 

Z,Sp,I,2,B 

0.7 

52 

2.0 

0.35+0.06 

21/24 

15.6+0.3 

21.8-22.5 

4 

4U 1758-25 

GX 5-1 

Z,Sp,I,-,B 

4.0 

73 

1.6 

0.73+0.05 

48/23 

12.9+0.2 


5 

4U 1758-20 

GX 9+1 

A,Sp,I,2,B 

1.2 

40 

1.7 

0.58+0.10 

35/23 

10.5+0.2 


6 

4U 1811-17 

GX 13+1 

A,Sp,I,-,B 

0.8 

23 

2.0 

0.63+0.06 

57/23 

11.6+0.2 

592.8 

7 

4U 1813-14 

GX 17+2 

Z,Sp,I,2,B 

2.3 

47 

1.4 

0.51+0.06 

31/24 

11.7+0.2 


8 

4U 1728-16 

GX 9+9 

A,Su,I,-,B 

0.3 

17 

1.5 

0.85+0.11 

35/23 

10.3+0.2 

4.2 

Type II 

9 

4U 1636-53 


A,Su,II,l,D 

0.2 

14 

4.0 

0.97+0.07 

35/24 

25.9+0.5 

3.8 

10 

4U 1735-44 


A,Su,II,-,D 

0.4 

12 

2.4 

0.85+0.07 

51/24 

22.4+0.4 

4.65 

11 

4U 1744-26 

GX 3+1 

A,Sp,I,-,B 

0.3 

21 

3.0 

1.13+0.10 

21/23 

22.8+0.5 


12 

4U 2142-b38 

Gyg X-2 

Z, 

1.4 

38 

3.7 

1.17+0.07 

67/25 

24.1+0.4 

236 

Type III 

13 

4U 1705-44 


A,Su,-,2,D 

0.13 

13 

30 

1.22+0.12 

76/24 

67.5+1.2 


14 

4U 1728-33 

GX 354-0 

A,Su,II,l,D 

0.1° 

6 

7.5 

1.79+0.11 

55/24 

41.7+0.8 


15 

4U 1820-30 


A,Su,II,-,- 

0.6 

19 

6.0 

1.59+0.11 

34/24 

32.4+0.6 

0.19 

Other sources 

16 

4U 1956-b35 

Gyg X-1 

BH 

0.7“ 

30 

10 

1.04+0.11 

129/25 

54.0+0.9 

134.4 


a Classification of the system in the various schemes (see text): Z, A=atoll, Sp=supercritical, Su=subcritical 
B=bulge, D=disc, I=class I, II=class II, l=one spectral component 2=two spectral components BH=black-hole source 
b From Christian & Swank (1997) in the energy range 0.5-20 keV 
c From Narita et al. (2001) in the 1-10 keV bandwidth 
d From Belloni et al. (1996) in the 1-30 keV bandwidth 
e Mean ASM count s“^ in the energy range 1.3-12.2 keV 

/am max/Imin 

g Best-fit power-law index of the power spectra. Errors are the Ict confidence interval 
h value and number of degrees of freedom of the fits to the power density spectra 
i Fractional (%) rms amplitude obtained from the light curves as the ratio of the variance over the mean 
j From Liu et al. (2001) 


the spectrum at low frequencies (below ~ 1 Hz) and a 
power law plus exponential cut off with index 0-0.8 (atoll 
sources), 0 (Z sources) and r'cut = 0.3 — 25 Hz (atoll 
sources), t'cut = 30 — 100 Hz (Z sources) at higher fre¬ 
quencies (van der Klis 1995 and references therein). The 
strength of these two components correlates with the posi¬ 
tion of the source in the colour-colour diagram. On top of 
this continuum several types of quasi-periodic oscillations 
are seen (van der Klis 1994a, 1994b; Wijnands & van der 
Klis 1999; Belloni et al. 2002). On long time scales (> few 
tens of days) the studies on LMXB have concentrated on 
the search for periodicities (Smale & Lochner 1992, Kong 
et al. 1998). In this work we investigate the long-term ape¬ 
riodic variability of 9 atoll and 6 Z persistent LMXBs by 
comparing the characteristics of their power density spec¬ 
tra in the frequency range 1x10“^ — 5xlO“®Hz and the 
patterns of variability of their light curves. 

2. Data analysis 

We have analysed the light curves obtained by the All Sky 
Monitor (ASM) on board the Rossi X-ray Timing Explorer 


(RXTE) of all persistent neutron-star systems showing an 
average count rate above 5 c/s (Table |]). The data were 
retrieved from the Definitive Products Database. The time 
span by the observations is about 5.7 years, from February 
1996 to October 2001 (JD 2,450,130-2,452,200). The ASM 
consists of three wide-angle (6° x 90°) shadow cameras 
(SSCl-3) equipped with position-sensitive Xenon propor¬ 
tional counters with a total collecting area of 90 cm^. The 
ASM scans ~ 80% of the sky every ^ 90 minutes in a se¬ 
ries of dwells of about 90 s each. Any given X-ray source 
is observed in about 5-10 dwells every day. The ASM is 
sensitive to X-rays in the energy band 1.3-12.2 keV. For 
more information on the ASM see Levine et al. (1996). 

In order to estimate the power density spectrum (PDS) 
of each source, we used 1-day binned light curves. Gaps 
due to detector failure or lack of data were filled by lin¬ 
ear interpolation, adding appropriate random noise. In 
order to avoid having a strong dependence on interpo¬ 
lation, which might introduce undesired features in the 
power spectrum, data interpolation was done only when 
the number of consecutive missing points was smaller than 
2% of the total number of points in the light curve. These 
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short gaps are randomly distributed over the whole light 
curve. The average number of missing points per gap was 
typically 2.5 for the 1-day rebinned light curves. 

The long gaps in the light curves, i.e. those occurring 
when no data points exist during an interval of time longer 
than 2% of the total length of the observation, divide the 
light curve into segments of different duration. In order 
to extract the maximum amount of information from the 
light curves we used the discrete (slow) Fourier transform 
on each observational segment of data. Since the data seg¬ 
ments had different lengths, the resulting power density 
spectra covered different frequency ranges. The PDS of 
the individual segments were merged together, sorted in 
order of increasing frequency, and rebinned in frequency 
to have at least 40 points per bin. The PDS were white 
noise subtracted and normalized such that their integral 
over a frequency range i^i — V 2 gives the squared fractional 
rms variability of the light curves due to variations on 
time scales from to (Belloni & Hasinger 1990, 
Miyamoto et al. 1991). 

3. Results 

Figure]^ shows the power density spectra (PDS) of all the 
sources studied in this work, together with the best-fit 
power-law model and residuals. The residuals represent 
the number of cr that the observations deviate from the 
model. For the sake of comparison we have also analysed 
the ASM light curve of the black-hole source Cyg X-1 in 
the same way (for details on the aperiodic variability of 
Cyg X-1 at very low frequencies see Reig et al. 2002). All 
sources show a red-noise dominated power density spec¬ 
trum (PDS). No flattening at low frequencies is seen in any 
of the PDS. Although a single power law did not give ac¬ 
ceptable fits in some cases, it provides a convenient way to 
compare the variability of the spectral parameters among 
the different types of systems. Characterising the LMXBs 
in terms of the shape of the power density spectrum alone 
may give rise to ambiguities since intrinsically different 
time series can produce similar power density spectra (e.g. 
Press 1978). That is to say, two sources may display the 
same PDS even though their variability mechanism may 
not be the same. Therefore we have also examined the X- 
ray variability of the sources in the time domain by obtain¬ 
ing the root mean square (rms) and the intensity distribu¬ 
tion function (IDF). The rms amplitude gives a measure 
of the variability of the source and was computed from the 
light curve as rms = jx^, where x is the mean count 
rate and cr^ = — tTg^p is the difference between the ob¬ 

served variance, = '^^{xi — x^/N, and the expected 
variance, Cg^p = '^i experimental er¬ 

rors, and N is the total number of points). Figure]^ shows 
the light curves of the sources analysed here. The IDFs 
(Fig. ^ show how often a given (normalised) intensity oc¬ 
curs. The abscissa represents the intensity of the source 
normalised to its mean value in steps of 0.1, whereas the 
ordinate gives the number of times (i.e. the frequency) 
that certain values of the normalised intensity occur. The 


IDF plots also give a measure of the source variability am¬ 
plitude. Narrow IDFs indicate that the source is roughly 
stable, showing low-amplitude variations, whereas highly 
variable sources will show many populated bins in their 
IDF plots, hence resulting in wide IDFs or narrow IDFs 
with extended tails. 

Figure ^ shows the rms of the light curves and 
the best-fit PDS power-law indices. By comparing these 
power spectral and timing parameters LMXBs can be 
divided into three groups. The first group includes all 
Z sources, excluding Cyg X-2, plus the atoll sources 
GX13-I-1, GX9-I-1 and GX9-I-9. These are bright sources 
showing low-amplitude and fast variations (on time scales 
of ^ 1-2 days). As a result, their rms values are small, 
their PDSs are flat and their IDFs appear to be nar¬ 
row and tend to display positive tails. The second group 
comprises the atoll sources 4U 1636-53, 4U 1735-44 and 
GX3-I-1 and the Z source Gyg X-2. These sources show 
similar ASM light curves as the previous group (i.e. low- 
amplitude variations on time scales of ~ 1 — 2 days) but 
with an extra long-term modulation. This modulation pro¬ 
duces more power at low frequencies giving rise to steeper 
power spectra and increasing their rms variability. Their 
IDFs are also broader when compared to the IDFs of the 
sources in the previous group and tend to exhibit negative 
tails. The third group contains those systems that cannot 
be placed in neither of the two previous groups, namely 
4U 1705-44, 4U 1820-30 and GX354. Although they do 
not form as homogeneous a group as Type I and II sources 
they display the highest amplitude variations, the steepest 
power spectra, and the broadest IDFs. Quasi-periodic os¬ 
cillations seem to be a characterising feature of this group. 
For the purpose of discussion we will refer to these groups 
as Type I, Type II and Type III, respectively. 

For comparison, all figures include data of the black- 
hole source Gyg X-1. Its PDS slope agrees with that of 
Type II objects, although it is also consistent with that 
of the Type III source 4U 1705-44. Its light curve and 
IDF are, however, rather peculiar (note the long tail above 
///mean = 1-5). 

Since sources belonging to the same group show very 
similar characteristics we have obtained a mean power 
density spectrum for each group by averaging the PDS of 
the individual sources pertaining to that group (Fig. ||). 
The PDS of Gyg X-1 has also been plotted (filled circles). 
Each group can then be distinguished on the basis of the 
slope of the PDS and the amplitude of variability. The 
power-law index a and the fractional amplitude of vari¬ 
ability rms calculated by integrating over the power-law 
model in the frequency range 1 x 10“^ — 5 x 10“® Hz of the 
mean PDS for each group of systems are: a = 0.62 ± 0.03, 
rTOS=11.5±0.3% for Type I sources; a = 1.09 ± 0.05, 
r?7is=15.0±0.4% for Type II sources and a = 1.61 ± 0.08, 
rTOS=29.5±0.6% for Type III sources. 

Type I and Type II systems exhibit similar PDS, ex¬ 
cept for the fact that Type II sources contain less (more) 
power at higher (lower) frequencies than the Type I sys¬ 
tems, which agrees with the characteristic shape of the 
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10 ^“ 10 ^^ 10 “’ 10 ^“ 10 ^^ 10 “'' 10 “" 10 “" 


Frequency (Hz) 


Fig. 1. Power density spectra of the sources studied in this work and the best-fit power law. The residuals represent 
the difference between the observational points and the model, expressed as number of cr. 


light curves, i.e. fast variations in Type I LMXB and long¬ 
term trends in Type II LMXB. Type III power spectra 
show the largest amplitude variations at all frequencies 
(except at the highest ones). This again is in agreement 
with the flaring like appearance of their light curves. 

Finally, we find a correlation between source luminos¬ 
ity, Lx, and rms, in the sense that more luminous systems 
tend to be less variable. In Fig.^, we plot the rms vs source 
luminosity. An anticorrelation between these two quanti¬ 
ties is apparent in this figure. This result is verified when 
we used the Kendall’s r nonparametric statistic (Press et 
al. 1992) in order to investigate, quantitatively, whether 
the apparent anti-correlation between the two quantities 
is significant or not. We find t = —0.56; the probability 
that we would get this value by chance if Lx and rms 
were uncorrelated is ~ 3.6 x 10“^. We fitted a straight 
line to the log(Lx) — log(rms) plot using the least-squares 
bisector line method (Isobe et al. 1990). The results show 
that rms oc L(^^, with (3 = 0.6 ± 0.1. 


4. Discussion 

We have investigated the aperiodic X-ray variability of 
all persistent neutron-star systems showing an average 
count rate above ASM RXTE 5 c/s over a period span¬ 
ning about 5.7 years. A new scheme is proposed in which 
LMXBs can be separated into three groups according to 
the shape of their power density spectra at very low fre¬ 
quencies and the pattern of the long-term X-ray variabil¬ 
ity of their light curves. As classifying parameters we have 
taken i) the fractional amplitude of variability rms from 
the light curves, ii) the slope a of the best-fit power law 
from the power spectra and, in a more qualitative way, Hi) 
the shape of the intensity distribution functions. 

Type I sources show lower rms and flatter power 
spectra {rms < 20%, a < 0.9) than Type II sources 
(20% < rms < 30%, 0.9 < a < 1.2), and these in turn, 
have lower rms and flatter power spectra than Type III 
{rms > 30%, a > 1.2). The intensity distribution func- 
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Time (days) 


Fig. 2. ASM light curves of the sources studied in this work. Each point represents 1 day. 

tions of Type I sources tend to show positive tails and Several schemes have been put forward for the classi- 
are narrower than those of Type II sources, which tend fication of LMXBs: Parsignault & Grindlay (1978) sepa- 
to show negative tails. In Type III systems the intensity rated LMXBs into class I and II depending on whether or 
distribution functions are even wider. not a direct correlation between the source temperature 

(obtained by htting an exponential to the energy spec¬ 
trum) and the source intensity exists. Ponman (1982) dis- 
It should be noted that the use of just one of these tinguished between subcritical and supercritical LMXBs 
parameters does not generally suffice to characterise un- depending on whether or not the X-ray luminosity of the 
equivocally a system. It is the combination of the spectral system is below or above a certain critical luminosity, ap- 
and timing parameters together with the extra support proximately half the Eddington limit. White & Mason 
from the characteristic shape of the intensity distribution (1985) classified LMXBs according to the number of corn- 
function which allows a more robust way of characterising ponents needed to fit their energy spectra, hence distin- 
LMXBs. For example, the slope of the PDS in GX9-I-9 is guishing between one-component (a single power law with 
similar to those found in Type II sources. However, the an exponential cutoff) and two-component (also including 
rms is about 10 sigmas smaller, more in accordance with a blackbody function) systems. Based on the age, the lo- 
Type I sources as its IDF is. The light curve exhibits a cation in the Galaxy and whether or not a low-energy 
long-term modulation similar to but lower in amplitude blackbody component was required to fit the energy spec- 
than those of Type II systems. Likewise, the power-law tra Naylor & Podsiadlowski (1993) divided the LMXBs 
index of Cyg X-2 is consistent with the Type III source into bulge and disc sources. According to the pattern that 
4U 1705-44 but its rms amplitude and shape of the IDF the source traces out in the colour-colour X-ray diagram 
resemble those of Type II sources. 
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Fig. 3. Distribution of the ASM intensity (normalised to the mean value) of the sources studied in this work. 


and the type of noise components in the 10“^ — 10^ Hz 
power density spectra, Hasinger & van der Klis (1989) cat¬ 
egorised LMXBs into Z and atoll sources. Note that the 
resulting number of different types is always two. 

In general, the agreement between these classification 
schemes is good (see Col. 4 in Table . One group would 
contain Parsignault & Grindlay’s class I, Penman’s super¬ 
critical, White & Mason’s two-component spectra, Naylor 
& Podsiadlowski’s bulge and Hasinger & van der Klis’ Z 
sources. The other group would include class H, suberiti- 
cal, one-component, disc, atoll sources. Nevertheless, each 
scheme has some controversial sources that detach from 
the norm. For example, GX349-I-2 (Kuulkers & van der 
Klis 1998) and GX13-I-1 (Homan et al. 1998) exhibit both 
Z and atoll characteristics. GX9-I-9 has a much shorter or¬ 
bital period than the rest of its group (bulge sources in the 
Naylor & Podsiadlowski’s scheme) making it very unlikely 
that this system contains a subgiant secondary as has 
been proposed for the bulge sources (Bandyopadhyay et 
al. 1999). GX340-I-0 did not easily fit in Ponman’s scheme 


(Ponman 1982), showing characteristics of both, suberiti- 
cal and supercritical sources. 

When the properties of the long-term variability are 
taken into account the classification of LMXBs based on 
the properties of the aperiodic variability at higher fre¬ 
quencies is mixed up. For example, Gyg X-2, which is 
considered as the prototype of Z sources, shows a much 
steeper power spectral continuum (a = 1.2) and higher 
amplitude of variability {rms=24%) than the rest of the 
Z sources, whose power-law index and rms distribute 
around values 0.7 and 13%, respectively. Likewise, atoll 
sources do not seem to form an homogeneous group ei¬ 
ther. Among the atoll sources the power-law index and 
rms variability vary over a wide range from 0.6 to 1.6 and 
10-60%, respectively. Some atoll sources such as GX13-I-1 
or GX9-I-1 are more similar to Z sources than to other 
members of the atoll group. 

It might very well be that the classification of LMXBs 
in terms of a discrete number of groups is an oversimplifi¬ 
cation and that LMXBs constitute an homogeneous group 
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of sources in which some unknown parameters change con¬ 
tinuously to give rise to their differences. 

One of these parameters could be the mass accretion 
rate. The increase of the X-ray luminosity from Type III 
to Type I objects may be interpreted as an increase in the 
accretion rate, probably due to the different nature of the 
secondary. It was already pointed out by Hasinger & van 
der Klis (1989) that the differences between Z and atoll 
sources may be due to the different luminosity class of the 
mass-donating stars, with Z sources having evolved sec¬ 
ondaries and atoll sources main-sequence stars. Infrared 
spectroscopic measurements (Bandyopadhyay et al. 1999) 
show evidence indicating that not only the luminosity 
class but also the spectral type might be different so that 
Z sources would contain evolved and earlier than G5 stars 
and atoll sources late-type (K or M), either evolved or 
main-sequence, stars. Since the mass-loss rate scales with 
luminosity and radius (Reimers 1975), it is natural to ex¬ 
pect that the mass accretion rate will be higher in systems 
which have an evolved star as an optical companion (like 
the Type I systems). 

If the X-ray variations are caused by variations of the 
accretion rate which propagate toward the inner regions of 
the accretion disc, then in Type I systems, the accretion 
disc should be able to transfer even the fastest of these 
variations (faster even than the diffusion time-scale at the 
outer part of the disc). This can be achieved if the accre¬ 
tion disc is very small in size (perhaps the inner part is 
missing due to the presence of a relatively strong magne¬ 
tosphere), or if it is geometrically thick (e.g. Churazov et 
al. 2001). 

As the accretion rate decreases, the disc becomes ge¬ 
ometrically thin, allowing only the slowest (and proba¬ 
bly largest amplitude) mass-rate variations to propagate 
toward the inner part. This could explain the very slow 
modulation in the light curves of Type II objects. Type 
III objects could simply be those Type II objects which 
have: a) the closest companion (note that the only value 
of the orbital period known for a Type III system is much 
shorter than any other system, see Table |^) and b) the 
weakest magnetic field. In this case, if the accretion disc is 
indeed smaller in size (due to the proximity of the compan¬ 
ion), then higher frequency oscillations could also propa¬ 
gate from the outer boundary toward the inner parts of 
the disc. At the same time, the larger amplitude mass 
accretion rate variations could penetrate the weak mag¬ 
netosphere and extend all the way to the surface of the 
neutron star, giving rise to the flaring appearance of their 
light curves. 

This scenario could also explain the dependence of rms 
on X-ray luminosity that we observe (Fig. ||). Suppose 
that the observed variations are indeed caused by varia¬ 
tions of the mass accretion rate. These variations could be 
thought of as ’’flares”, which originate at the outermost 
part of the disc and propagate toward the innermost re¬ 
gions. Suppose now that the shape and size of the flares 
are fixed in all systems, and that only the number of flares 
per unit time, say N, differ among the various systems. In 



2 4 6 8 10 12 14 16 

source number 


Fig. 4. rms amplitude of the light curves and power-law 
index of the power spectra of the individual sources (see 
Table 1^. Type I sources are represented by filled squares. 
Type II sources by open squares and Type III sources by 
stars. The filled circle corresponds to Cyg X-1. The errors 
in rms amplitude are smaller than the symbols. 


this case, assuming there is not a significant constant, un¬ 
derlying component. Lx oc N, while rms oc N/L'^, and 

so rms oc inconsistent with the observed relation 

— 0 6 ^ 

rms oc ' . Therefore, both N and the typical flare size 
scale (i.e. duration) should change with luminosity. If, as 
we move from Type I to Type II and Type III objects, only 
the longest amplitude variations propagate in the disc (i.e. 
N decreases) and these variations last longer (i.e. the flare 
size scale increases), then a relationship between rms and 
Lx, similar to the one shown in Fig. || could be expected. 

It is worth making a comparison between Cyg X-1 
and neutron-star systems. The power spectra of Type I 
sources are significantly different from that of Cyg X-1 
both in slope and variability. Type II sources, although 
they exhibit similar noise shape, are less variable than 
Cyg X-1. In Type III sources the rms is comparable to 
that of Cyg X-1 but have a steeper slope in the PDS. 
The differences between the properties of Cyg X-1 and 
some of the LMXB such as 4U 1705-44 (whose power- 
law index and rms amplitude are similar to those of Cyg 
X-1) are smaller than between some LMXBs themselves. 
We conclude that the physical processes that cause the 
low-frequency X-ray variations in X-ray binaries are de¬ 
termined by the conditions of matter in the surroundings 
of the compact object rather than by its nature (neutron 
star or black hole). 
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Fig. 5. Comparison of the power spectra of low-mass X- 
ray binaries with that of Cyg X-1. The power spectra are 
averaged spectra of the individual sources that belong to 
each group. Errors are absorbed by the size of the points. 
Symbols are as in Fig. |[ 
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Fig. 6. Fractional amplitude of variability {rms) as a func¬ 
tion of X-ray luminosity in the energy range 0.5-20 keV, 
except for the points marked with an arrow, GX354 and 
Cyg X-1, for which 1-10 keV and 1-30 keV bandwidths 
apply, respectively. The dashed line represents the best fit 
to the data using the least-squares bisector line method. 
Symbols are as in Fig. |[ 

5. Conclusions 

Based on the properties of the long-term aperiodic vari¬ 
ability we have divided low-mass X-ray binaries into three 
categories. Systems in each group distinguish themselves 
by the shape and amplitude of variability of their power 
density spectra, light curves and intensity distribution 
functions. It seems plausible that Type I systems har¬ 
bour evolved and possibly early type companions, whereas 
Type II and III systems contain later type main-sequence 
stars. The mass loss of the companions increases from 


Type III to Type II to Type I systems, hence the corre¬ 
sponding increase in luminosity. Differences in the accre¬ 
tion disc structure (thin vs thick) could explain the lack of 
fast variations in Type II systems. Type III systems could 
also be affected by the weak magnetic field, and perhaps 
the small size of the disc. Our results imply that, after all, 
it is not just the accretion physics that dictates the ob¬ 
served characteristics of the systems. Other factors, like 
the size of the disc and/or the presence of weak magnetic 
fields, probably affect the low frequency X-ray variations 
in these systems. Although the separation of LMXBs into 
these three groups is useful for the purpose of investigating 
their general aperiodic properties, it is likely that LMXBs 
constitute an homogeneous group of sources, whose prop¬ 
erties vary in a continuous, rather than a discrete, way. 
Whether the differences in the aperiodic variability be¬ 
tween Cyg X-1 and LMXBs are due to the different na¬ 
ture of the primary star - a B supergiant in opposition to 
late G-K type stars - or intrinsic variations in the vicinity 
of the compact star is unknown. It would be interesting to 
investigate the long-term variability of black-hole systems 
with late-type companions. However, owing to the tran¬ 
sient nature of these systems such studies are very difficult 
to perform. 
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